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ABSTRACT 

The recursion relations that were proposed in [2] for implementing vector extrapolation 
methods are used for devising generalizations of the power method for linear operators. 
These generalizations are shown to produce approximations to largest eigenvalues of a linear 
operator under certain conditions. They are similar in form to the quotient-difference algo- 
rithm and share similar convergence properties with the latter. These convergence properties 
resemble also those obtained for the basic LR and QR algorithms. Finally, it is shown that the 
convergence rate produced by one of these generalizations is twice as fast for normal opera- 
tors as it is for non-normal operators. 


•Work funded under Space Act Agreement C99066G. 
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Let B be a normed linear space over the field of complex numbers, and denote the norm associ- 
ated with B by 1 1 • 1 1 . In case B is also an inner product space, we adopt the following convention for 
the homogeniety property of the inner product: for y ,z € B and a,p complex numbers, the inner pro- 
duct (v) satisfies (ay ,(3z) = aP(y ,z). The norm in this case is the one induced by the inner product, 
i.e., ifxeB, I lx II = ^(x,x). 

Let x, , i =0,1,..., be a sequence in B. We shall assume that 


asm->«. 

i=i 


( 1 . 1 ) 


Here v 1( v 2 , . . . , are linearly independent vectors in B. X, are distinct scalars ordered such that 


IXjI > IXjI > IX3I >-, 


( 1 . 2 ) 


and satisfy 


X, * 0, i = 1,2,...; X, * Xj if 1 *j . 


(1.3) 


Also there can be only a finite number of X, ’s having the same modulus. The meaning of (1.1) is that 
for any positive integer N there exist a positive constant K and a positive integer m 0 that depend only 
on N, such that for every m>m 0 , 

Ux m - xV-Xril <*lX w r. (1.4) 

1=1 

Sequences of the form described above arise naturally, for example, in the matrix iterative pro- 


cess 


Xy +1 =AX ; , j =0,1,..., (1-5) 

in which Xq is an arbitrary column vector, and A is a square matrix whose nonzero eigenvalues have 
only associated eigenvectors and no principal vectors. In this case X^Xj, . . . , are some or all of the 
nonzero eigenvalues of A , depending on the spectral decomposition of Xq, and vj,V 2 , . . . , are propor- 

00 

tional to Its corresponding eigenvectors. Also the infinite sum in (1.1) is now finite, and the 


asymptotic relation there becomes an equality. 
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In the next section we present three recursive techniques, which we designate the qd-MMPE, 
qd-TEA, and qd-MPE algorithms, that are based solely on the vector sequence *o**i»*2» • • • . and that 
can be used for obtaining approximations to X^XiXi , .... in this order. These techniques generalize 
the well known power method and resemble the quotient-difference (qd) algorithm. In Section 3 we 
analyze the convergence properties of these techniques and show that they indeed behave like the qd 
algorithm. We recall that the qd algorithm is used in approximating the poles of meromorphic func- 
tions in general, and the zeros of polynomials in particular. For detailed descriptions of the qd algo- 
rithm, see, for example, Henrici [3, Chap. 7] and Householder [4, Chap. 3]. The convergence rates 
derived for the methods of the present work are also very similar to those derived for the basic LR 
and QR algorithms, see, for example, Parlett [5]. 

In Section 4 we will analyze one of the recursive techniques of Section 2, namely the qd-MPE 
algorithm, in conjunction with vector sequences x m ,m =0,1,..., that are generated by iterating with a 
normal operator, and we will show that the rate of convergence for this case is twice that obtained for 
an arbitrary non-normal operator that has the same spectrum. 

Finally, we mention that the recursive techniques developed in this work are based on the recur- 
sive algorithms that were developed by Ford and Sidi [2] for implementing some vector extrapolation 
methods, such as the modified minimal polynomial extrapolation (MMPE), the topological epsilon 
algorithm (TEA), and the minimal polynomial extrapolation (MPE), from which their names are 
derived. The methods of proof relating to the convergence of these techniques are similar to those that 
were used in Sidi, Ford, and Smith [8] and in Sidi [6] in the analysis of the above mentioned vector 
extrapolation methods, and in Sidi [7] in the analysis of some other recent extensions of the power 
method as they are applied to normal operators. We also mention that the scalar and vector epsilon 
algorithms of Wynn [9,10] can also be used as generalizations of the power method and produce 
results similar to the ones given in Section 3 of the present work. This has been shown by Brezinski 
[ 1 ]. 
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2. DEVELOPMENT OF THE ALGORITHMS 

In this section we follow very closely the developments and notation of [2]. 

In all three methods that we develop in this section we assume that the vector sequence 
x 0 JiJ 2 - • • • » is given. We also assume that a sequence of scalar quantities 1 l”, m,n >0, is given. 
These quantities are determined in different forms from the vectors x t . Their exact definition will be 
given later in this section. In addition, we denote by b an arbitrary sequence b 0 ,bi,b 2 ,— . Thus p/ 
stands for the sequence p<{,p./,p .£... . 



and/* ^(b) to be the determinant 

b n b/t+l ' ' ' t>n+k 

iC iCl ' ' ' iC* 

| ’ 
^ + +*-l ... ^ + +*- 

We also define 


Go 1 ” = 1 and f n 0 /n (b) = b n . 

Finally, we let 

T* ^{b ) = b )/G; +1 ' m . 

Note that by (2.3) and (2.4) 


( 2 . 2 ) 

(2.3) 

(2.4) 


rr(b) = b n . 


(2.5) 
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2.1 The qd-MMPE Algorithm 

Let 2 1,2 2 , . . . , be a sequence of linearly independent bounded linear functionals on B, and let 

K = Qm+ 1(*»). m > n - 0- (2-6) 

We recall that in case B is a complete inner product space each Qj has a unique representer q t in B in 
the sense that Qj(z) = (qj ,z ) for all z in B. 

As is shown in [2, Theorem 2.1], the quantities T£(b) = T£-°(b) satisfy the 3-term recursion rela- 
tion 

mb ) = m x (b ) - dtmi & ), ( 2 . 7 ) 

where 

= n_x 01 k ~ x ymi (n* _1 ) . ( 2 . 8 ) 

Note that d\ = from (2.5) and (2.8). Note also that T£(yJ) = 0 for 0 <j < k- 1. 

As will be shown in Section 3, for the vector sequences x 0 qc x jC 2 ’ • • - « described in the begin- 
ning of Section 1, 1 ld k — ^ ^ as n — ) °° under certain conditions. Thus the developments above pro- 
vide us with a technique that can produce approximations to X^Xj, . . . , in this order. This technique 
is summarized below: 

qd-MMPE Algorithm ( 2 -9) 

Given the sequence Xo>*i>* 2 «— . compute |i™, m ,n >0, by (2.6), and set Tq(\J) — For 
k = 1,2,..., use (2.8) to compute d k , and then use (2.7) to compute 7*(jj J),j > k . 

We note that d k is constructed from the £+1 vectors x n oc n+l ,...oc n+k . 

Finally, we observe that \!d\ = p n V^ 0 = Q \(x n+x )lQ i(*„). which is exactly the approximation 
provided by the power method for Xj, the largest eigenvalue in modulus of a linear operator, provided 
this eigenvalue is simple, and satisfies I Xj I > 1 X^ I . In this sense the qd-MMPE algorithm general- 
izes the power method. 
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2.2 The qd-TEA Algorithm 

Let Q be a bounded linear functional on B, and let 

iC = 2 (*«+„). m,n>0. (2.10) 

Again, in case B is a complete inner product space Q has a unique representer q in B in the sense that 
Q(z) = ( q ,z) for all z in B. The rest is exactly as in the qd-MMPE algorithm. We only have to observe 
that d k now is determined by the 2k vectors x H JC n +i>*—x n + 2 ji-i- We also observe that for k = 1 the qd- 
TEA and qd-MMPE algorithms are the same when Q=Q\. Consequently, as with the qd-MMPE 
algorithm, also with the qd-TEA algorithm the power method is obtained for £=1. Hence the qd-TEA 
algorithm also generalizes the power method. 


2J The qd-MPE Algorithm 

The scalar quantities j are now defined by 


1 l? = (x m qc„), m,n> 0. (2.11) 

As is shown in [2, Theorem 3.1], the quantities T n k {b) = T n k *{b) and T n k {b) = T^~\b) satisfy 
the coupled recursion relations 


T n k (b) = n_ x (b)-dZf£l{b), n> 0, 
f k (b) = T k _i (b) - (1 b ), n > 1, 


( 2 . 12 ) 


where 


d n k = r;_! (vr+t-'yffi (p' ,+A:_1 ), n > o, 

d n k = TU 0* " _1 )/2*_ + i 1 Ot"- 1 ), n > 1. 


(2.13) 


In fact, there exists a very interesting 4-term (lozange) recursion relation among T k+l (b), T k (b), 
T k +l (b), and T k + X ( b ), see [2]. Note that d\ = p"/p, B + i, n > 0, and d" =p n "" 1 /|i.«+i 1 . n > 1, from (2.5) 
and (2.13). Note also that T k (\i j ) = 0 for n<j < n+k- 1, and f k {yJ) = 0for«-l <j <n +k- 2. 

As will be shown in Section 3, for the vector sequences x 0 ^ijc 2 . described in the begin- 

ning of Section 1, 1 !d k -» \ k and 1 /d k -> X* a sn under certain conditions. Thus, the develop- 
ments above provide us with yet another technique that can produce approximations to Xi,X 2 , ...» in 
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this order. This technique is summarized below: 

qd-MPE Algorithm (2.14) 

Given the sequence x 0 jc 1 jc 2 ’ • • • » compute m,« > 0, by (2.11), and set 

Tq (p/) - \il = Tq( p/). For k = 1,2,..., use (2.13) to compute d* and and then use (2.12) 
to compute 7*(p. y )» j - n—l,j ^ n+k, and T£(pj), j < n-2,j > n+k—l. 

Observe that d£ and d" are constructed from the k+2 vectors x n _ijc n jc n +i,...jc n+k only. 

Finally, we observe that lid" = p" +1 /p,? = (x„ pc nJrX )l(x n ,x n ), which is exactly the Rayleigh quo- 
tient that provides an approximation to A.}, the largest eigenvalue in modulus of a linear operator, pro- 
vided this eigenvalue is simple, and satisfies 1^ I > I I . In this sense the qd-MPE algorithm gen- 
eralizes the power method (the Rayleigh quotient). It is known that in case the vector sequence 
• • • . is generated by a normal operator, the Rayleigh quotient converges to Xi twice as 
quickly as it does for a non-normal operator. 

So far we have shown that the three algorithms above generalize the power method for a linear 
operator. We would now like to explain why they are of the quotient-difference type. First, the algo- 
rithms are somewhat similar in form to the quotient-difference method. Second, as is shown in the 
next section (see (3.20) and (3.24)), the d£ are all expressed in terms of four determinants the way the 
quantities in the quotient-difference method are. Furthermore, again as will be shown in the next sec- 
tion, these determinants have asymptotic expansions very similar in form to those obtained for the 
determinants involved in the quotient-difference method. We note in passing that the determinants that 
appear in the quotient-difference algorithm are Hankel determinants. The same holds true for the qd- 
TEA algorithm as will be clear in the next section, although the indices of these determinants are dif- 


ferent in the two methods. 
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3. CONVERGENCE ANALYSIS OF THE ALGORITHMS 

We now state the main convergence results for the qd-type algorithms that were devised in the 
previous section. 

The following assumption is common to all three methods: 

I 1 > IX* I > IX* +1 I for some integer k > 1; X^®*. (3.1) 


Let us also define 



x k 


^•*+1 




J 

X 

s 


e* = max 

Thus, e* < 1. 

Theorem 3.1: In the qd-MMPE algorithm denote 

2i(vi) Qi(vj) 


, k> 1. 


S/ = 


G/(v,) Qjiyj) 


, j = 1,2,. 


Provided (3.1) holds, and 


we have 


S*_j ^ 0 and S* ^ 0 , 


l/d k = X k +0 (e*) as n ->». 
Theorem 3.2: In the qd-TEA algorithm provided (3.1) holds and 


(3.2) 


(3.3) 


(3.4) 


(3.5) 


n£(v,)*o, 

i=l 


(3.5) is satisfied. 


(3.6) 


Theorem 3 J: In the qd-MPE algorithm provided (3. 1) holds, (3.5) is satisfied. (3.5) is satisfied also 

when d k is replaced by d k . 

Note: In relation to Theorem 3.3 we realize that if the sequence x 0 ^i^ 2 . • • • - is generated by the 

matrix A in the way x ;+1 —Ax } , j = 0,1,..., x 0 being arbitrary, then the qd-MPE algorithm performs 
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asymptotically like the qd-MMPE and qd-TEA algorithms, although it does not require additional 
conditions like (3.4) and (3.6). 

In the proofs of Theorems 3. 1-3.3 we shall make use of the following results, which are of 
interest in themselves. 


Lemma 3.4: Let <Ji,a 2 , .... and £ lt £ 2 , ■ ... be two sequences of nonzero scalars, such that the o,- 

are distinct, and 


ICiI>IC 2 I* • (3.7) 

Assume, furthermore, that there can be only a finite number of ’s having the same modulus. Con- 
sider the determinant 


in which 


W = 


« t? 


■P-i 


II ^ i/ (^) ... u (*) . 

“r f 0 w r, 1 “r/- 1 


(3.8) 


zjjj being some scalars. Define 


u p,q ~ Z z 7>°/C/ as n °°t 

7=1 


(3.9) 


Z • 


z jui z ji2 ' ' ’ z j\s 
Z h.l Z ji2 Z jiS 


Z jr. 1 Z 7,.2 

Then y", for n -4 =», has the asymptotic expansion 




(3.10) 


V" - 


Z 


m 

p m i 


where V(a 1( . . . , a r ) is the Vandermonde determinant defined by 


(3.11) 
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V(a x a r ) = 


1 a x 
1 a 2 


• «r' 

• a r 2~ l 


*r- 1 


(3.12) 


1 a r ••• ay 

(When ~ is replaced by = in (3.9), - is replaced by = in (3.11) too.) If, in addition, 

I Cr I > I Cr+1 * 

then 



f ■" 

r 

■ f 

Ul 

*y 

1^,...^ V(Oi,<J2, . . . ,a r ) 

re; 

U-1 J 

l+O 

_ V 

Cr 

/ . 


as n- 


(3.13) 


(3.14) 


Lemma 3.4 can be obtained from [7, Lemmas 2.2 and 2.3], whose technique is a generalization 
of that employed in [8, Theorem 3.2]. 

Lemma 3.5: Let ^,£ 2 , .... be a sequence of nonzero scalars exactly as described in Lemma 3.4. 

Consider the determinant \j/” as given in (3.8), but with 






(3.15) 


i=l Jm 1 


zij being some scalars. Define 




Z l't/l Z ljl 


Z iij\ Z ixjj 


* Z, 


‘l/r 




Z *rj\ Z irjl ' # ' 2 ir]r 

Then for n has the asymptotic expansion 


(3.16) 


W~ 


711,...^ 




/>=! J 


Lp-i J 


(3.17) 


(When - is replaced by = in (3.15), - is replaced by = in (3.17) too.) If, in addition. 
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iCr I > lC-+il and Z};z'*0, 


then 


iv(Ci Cr)' : 


r 

2 n 


Cr +1 


nc, 

y=i 


1+0 

C, 

J . 


as n 


(3.18) 


( 3 . 19 ) 


Lemma 3.5 can be proved exactly as in [6, Lemma 3.2 and Theorem 3.1]. 

Proof of Theorem 3.1: Combining (2.8) with (2.4), (2.2), and (2.1), and recalling that 

T£(b ) = T£ ,0 (b), we obtain 




G n ,0 2,0 

k W-l 

G n+ 1,0 

k ^k - i 


(3.20) 


Now by (2.6) and (1.1) we have 


K~ Z<2m+i(V;)V as n 
y=i 


( 3 . 21 ) 


Substituting (3.21) in the determinant expression for G* ,0 , we see that G* ,0 is simply vy* of Lemma 
3.4 with - ji£^, thus z j<p =Q p {Vj), Cy = X/, and a } =X } . Furthermore, S ; in (3.3) is Z IA ^ in 
Lemma 3.4. Therefore, Lemma 3.4 applies, and we have 

" **♦! 


= V(k lt ...,X k ) 


m 

M . 


1 + 0 




as n —* 


( 3 . 22 ) 


Treating G£ +I -°,G£! i 1,Q , and G£!?'° similarly, and combining the results in (3.20), we obtain (3.5). □ 


Proof of Theorem 3 2: We observe that by (2.10) and (2. 1) 


K -'Ll Q (v,)W as n -+ =o. ( 3 . 23 ) 

f=i 

The rest of the proof is exactly the same as that of Theorem 3.1, with Q m+ \(Vj ) in the latter being 


replaced by Q (v.-)Xj" now. On account of this, Zj j now becomes 


TO(vy) 

i=i 


V(Xj X,). We 


leave out the rest of the details. □ 

Here we note that G £ ,0 for the qd-TEA algorithm is a Hankel determinant. Thus the qd-TEA 
algorithm resembles very strongly the qd algorithm. Furthermore, when the sequence x 0 jc yjci , .... 
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is one of scalars, then the qd-TEA algorithm can still be employed by picking the functional Q to be 
the identity operator, in which case v£=x m+n . Thus for scalar sequences satisfying all the conditions 
mentioned in Section 1 of this work the qd-TEA algorithm provides another form of the qd algorithm. 


Proof of Theorem 33: Combining (2.13) with (2.4), (2.2), and (2.1), and recalling that 

T£(b ) = 7X - n (b ) and f£(b ) = T^~\b ), we obtain 


Gr Gtf* 


dl = 


GS* 




Qjt+lfi-l Q* + l* 


(3.24) 


Substituting (2.11) in (2.1), and comparing with (3.8), we see that G i "’' ,+, = Vt with 
uj$ = = (Xn+s+p-i^)- Invoking now (1.1), we see that uj$ is of the form (3.15), with 

z«y = ( v,-,v,-)Xf and = X,. Also in Lemma 3.5 now becomes 



(Vi.Vj) 

••• (Vi.v*) 



f’ 

k _ 


. 

- 


r > 

Jfc _ 

7 W_ 

nv 


. 

. 

= 

nv 


i-1 , 


. 

. 




(v*,Vl) 

••• ( v*,v* ) 



(3.25) 


Since the vectors v lp v 2 , . . . , in (1.1) are assumed to be linearly independent, the Gram determinant 
U x of the vectors v lt ...,v t is positive. Thus 0. Consequently, Lemma 3.5 applies, and we 


obtain 


G£' n+ * = 


nv 

1-1 


X*)l : 


k 

In 

■ 

Xjt+i 

n' 

nv 

>-i 


1+6? 

X k 

J . 


as n 


(3.26) 


Invoking (3.26) in (3.24), we obtain the required result. □ 
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4. THE qd-MPE ALGORITHM FOR NORMAL OPERATORS 

As we have shown in Theorem 3.3, the convergence rate of the qd-MPE algorithm is the same as 
those of the qd-MMPE and qd-TEA algorithms, in general. When the vectors v i,v 2 , . . . , in the expan- 
sion (1.1) are orthogonal with respect to the inner product (•,•). however, the convergence rate of the 
qd-MPE method becomes twice that given in Theorem 3.3, as will be shown shortly in Theorem 4.1. 
When the vector sequence •••» * s generated by the matrix iterative process Xj +l =Axj, 
j = 0,1,..., with x Q arbitrary, then, as we mentioned in the introduction to this work, x m satisfies (1.1), 
with X, there being nonzero eigenvalues of A that possess only eigenvectors and no principal vectors, 
and v, there being proportional to the corresponding eigenvectors. Now when A is a normal operator 
it has only eigenvectors, and these form an orthogonal set with respect to the inner product. We recall 
that heimitian and antihermitian operators are normal. Thus Theorem 4.1 below applies to vector 
sequences generated by matrix iterative processes with normal matrices. 


Theorem 4.1: In the qd-MPE algorithm, provided (3.1) holds and the vectors v lf v 2 , .... in (1.1) 

satisfy 

(v; ,Vj ) = Zj 5 ij , (thus Zj >0), (4.1) 

we have 


lld k = X k +0 (e*") as (4.2) 

When Xj > 0 ,j = 1,2,..., (4.2) can be refined considerably. (In this case the condition in (3.1) is 
automatically satisfied.) If X k+ \/X k > X k /X k _i, then 

2a 


Md n k -X k -$ k 


'•jfc+l 


as n — » some |3* > 0, 


(4.3) 


i.e., the sequence {l/d k }^ converges to X k monotonically from below, and if X k+i /X k < X k /X k _\, then 

2a 


1 'dk~X k +Y k 


^*-1 


as n -> «, some y k > 0, 


(4.4) 


i.e., the sequence {l/d?}”^ converges to X k monotonically from above. P* and y k are given in the 
proof below. 
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AU the results above hold when d k is replaced by d k , though with different P* and y*. 


Proof: We saw in the proof of Theorem 3.3 that G£ A+S = y* with u$ = (x n+s+p -ijc n+<l ). Invoking 

now (1.1) and (4.1), we see that ufy is of the form (3.9), Zjj, = ZjXj +p ~ 1 , a y = Xj, and £ y - = I \j I 2 . 


Also Zj x in Lemma 3.4 becomes 


Zju-Jk 


Lp=i J 


m 

L /’ =1 J 


V(Xj X Jt ). 


(4.5) 


Substituting all this in (3.1 1), we obtain for n — » 00 


G2 A+S ~ £ 


Eft 

P - 1 


n * 

p = 1 . 


1 V(X h X jt )\ 


UK 

p = 1 


2n 


(4.6) 


Consequently, all the conditions of Lemma 3.4 are satisfied, and 

2*1 


G k A+s = 


nv 

/-! . 


rT z ; 


T&j 

7=1 


1+0 


v Jk+l 


2/i' 


as /z (4.7) 


I V(X x A*)l‘ 

Invoking (4.7) in (3.24), (4.2) follows. 

Now when the X } ’s are all real, from (4.6) and (4.1), we see that G£ A+S , s = 0,±2,±4,..., are all 
real and positive. When the Xj's are all positive G k J are positive for all s. In case all Xj's are posi- 
tive, (1.2) and (3.1) become 

> X 2 ^**‘^^—1 ^ ^ ■ (4.8) 

so that (3.1) is automatically satisfied now. The two most dominant terms in the summation in (4.6) 
are those with the indices = 1,2,...,/: and j\,...J k - \2,...,k-\,k+\. Thus 

as n -+ <*>, (4.9) 


G£ A+ * = 

f ■> 

k 

n z ; 

W(X U .. 

.,X k )\ 2 

k 

n^ +s 

1 +a* 

^*+1 

X, 

2n+s 

+ 0 

\ x M 
x k 

2/iV 


u-» J 



j mi 

- 

,x k 

s* J 


K 

V J 

J . 


where 


a k = 


z k + 1 

1 

>> 


z k 

/— N 

>* 

* 

I 

/ •, 

<< 

7 

<< 


z k+ l *-l 
z k j. 1 


X k +r-Xj 


X k ~Xj j 


>0. 


(4.10) 


Consequently, for n 
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2*1 


d?~ 


^ * 


1-KXt 


1+Ott 


'•i+l 


l+CCj 




2/i+n 


Vt+i 


l+ai._ 


*-i 




^■*-1 


2/i+Zl 


2n+f| 


1 

K 


2*1 


2 /i+n 


1+OCi. 


H+l 


H+l 




l- 


'•t-i 


**- 


j_ 

h 


2a 


1*HXju 


1- 


H+l 


H+l 


-a*-i 




A. 

^-*-1 


2*+n 


The rest of the proof can now be completed easily. We only note that P* = a k 


1- 


St+l 


7t=<X*-l 


1- 


'-*-1 


K- 


in (4.4). □ 


■k - 1 


(4.11) 


in (4.3), and 


-16- 


REFERENCES 


[1] C. Brezinski, Computation of the eigenelements of a matrix by the e-algorithm, Lin. Alg. 
Appl., 11 (1975), pp. 7-20. 

[2] W.F. Ford and A. Sidi, Recursive algorithms for vector extrapolation methods, Appl. Numer. 
Math., 4 (1988), pp. 477-489. 

[3] P. Henrici, Applied and Computational Complex Analysis, Vol. 1, John Wiley & Sons, New 
York, 1974. 

[4] A.S. Householder, The Numerical Treatment of a Single Nonlinear Equation, McGraw-Hill, 
New York, 1970. 

[5] B.N. Parlett, The LU and QR algorithms, in Mathematical Methods for Digital Computers, 
(A. Ralston and H.S. Wilf, eds.), Vol. II, John Wiley & Sons, New York, 1967, pp. 116-130. 

[6] A. Sidi, Convergence and stability properties of minimal polynomial and reduced rank extra- 
polation algorithms, SIAM J. Numer. Anal., 23 (1986), pp. 197-209. 

[7] A. Sidi, On extensions of the power method for normal operators, Lin. Alg. Appl., 
120 (1989), pp. 207-224. 

[8] A. Sidi, W.F. Ford, and D.A. Smith, Acceleration of convergence of vector sequences, SIAM 
J. Numer. Anal., 23 (1986), pp. 178-196. 

[9] P. Wynn, On a device for computing the e m (S H ) transformation, MTAC, 10 (1956), pp. 91-96. 

[10] P. Wynn, Acceleration techniques for iterated vector and matrix problems, Math. Comp., 16 
(1962), pp. 301-322. 



NASA 

National Aeronautics and 
Space Administration 


Report Documentation Page 


1. Report No. 


NASA TM- 102361 
ICOMP-89-25 


2. Government Accession No. 


4. Title and Subtitle 


Quotient-Difference Type Generalizations of the 
Power Method and Their Analysis 


3. Recipient’s Catalog No. 


5. Report Date 
September 1989 

6. Performing Organization Code 


7. Author(s) 

Avram Sidi and William F. Ford 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


8. Performing Organization Report No. 

E-5077 


10. Work Unit No. 


505-62-21 


1 1 . Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


[14. Sponsoring Agency Code 


15. Supplementary Notes 

Avram Sidi, Institute for Computational Mechanics in Propulsion, Lewis Research Center (work funded under 
Space Act Agreement C99066G) and Technion— Israel Institute of Technology, Computer Science Department, 
Haifa 32000, Israel; William F. Ford, NASA Lewis Research Center. Space Act Monitor: Louis A. Povinelli. 

16. Abstract 

The recursion relations that were proposed in 4 4 Recursive Algorithms for Vector Extrapolation Methods” by 
W.F. Ford and A. Sidi (Appl. Numer. Math, 4 (1988), pp. 477-489) for implementing vector extrapolation 
methods are used for devising generalizations of the power method for linear operators. These generalizations are 
shown to produce approximations to largest eigenvalues of a linear operator under certain conditions. They are 
similar in form to the quotient-difference algorithm and share similar convergence properties with the latter. These 
convergence properties resemble also those obtained for the basic LR and QR algorithms. Finally, it is shown 
that the convergence rate produced by one of these generalizations is twice as fast for normal operators as it is for 
nonnormal operators. 


17. Key Words (Suggested by Author(s)) 

Eigenvalue problem 
Power method 
Extrapolation methods 
Quotient-difference algorithm 


19. Security Classif. (of this report) 
Unclassified 


18. Distribution Statement 

Unclassified — Unlimited 
Subject Category 64 


20. Security Classif. (of this page) 
Unclassified 


21. No of pages 
18 


22. Price* 


NASA FORM 1626 OCT 86 


For safe by the National Technical Information Service, Springfield, Virginia 22161 






